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The equation for gamma-ray degradation spectra is reduced to a form suitable for a 
simple approximation. The approximation is valid in a wavelength range between two 
ComptoD wavelengths larger than the wavelength of the source and the region where photo- 
electric absorption becomes very intense; this range is wide for low 7* materials. Input data 
are provided for application to various materials and source energies. The results of sample 
calculations are shown and compared with those obtained by standard numerical calculation. 



In the course of gamma-ray penetration, the spec- 
trum of multiply scattered, degraded, secondaries 
approaches a shape that depends only on the prop- 
erties of the medium. 1 This limiting spectrum is 
that which arises when a gamma-ray source is 
distributed uniformly throughout a medium and the 
only relevant factor is the accumulation at each point 
of secondaries of various energies. When degrada- 
tion spectra were first calculated numerically, 2 M. H. 
Johnson pointed out, in a private discussion, that 
the main features of the degradation spectrum at low 
energies are determined by very simple arguments. 
A formalism of degradation theory, starling from 
this remark and from the corresponding neutron 
theory, 3 was developed later. 4 But this theory was 
not applied in actual gamma-ray calculations, 
because these were being carried out routinely by 
direct numerical computation. 5 Lately, the exten- 
sion of calculations to degraded photon energies ^20 
kev has made it desirable to try to evaluate the simple 
approximate analytical formulas and compare the 
results with those obtained by numerical integration 
of the basic equations. These low photon energies 
are of interest, e. g\, in LiH. This paper contains a 
description of the procedure, beginning with the 
qualitative analysis, together with a discussion of the 
results of the comparison. 

If successive steps of degradation by Compton 
scattering, characterized by wavelength shifts 
£=X— X', have approximately equal statistical dis- 
tribution (see footnote 4) with a mean (£), the prob- 
ability for each photon to attain a wavelength be- 
tween X and X + AX at the end of one of the steps of 
degradation is AX/(£). 6 If the photon does "land" 
in this interval, it w r ill "spend" there a mean free 
path l/AWt(X), where ju 8catt is the probability of 
scattering per unit pathlength and coincides with the 
total /z(X) when absorption is disregarded. For a 
source strength of S photons per unit volume and 



I i See/c. g., U. Fano, Nucleonics 11, Nos. 8 and 9, 55 (1953). 

2 P. R. Karr and J. C. Lamkin, Phys. Rev. 76, 1843 (1949). 
n 3 See, e. g., R. Marshak, Rev. Mod. Phvs. 19, 185 (1947), and in particular, 
G. Placzek, Phvs. Rev. 69, 423 (1946). 

* U. Fano, Phys. Rev. 92, 328 (1963). 

5 Calculations of penetration by the moment method, as given by H. Goldstein 
and J. E. Wilkins, Jr., NYO Report 3075 (1954), include, under the name of 
"zero-tta moment", the spectrum discussed here. 

6 This problem is similar to the determination of the probability that a pedes- 
trian will step on a line drawn across a sidewalk. The probability equals the 
length of the pedestrian's foot divided by his mean step length. 



unit time, we thus expect at each point a flux of 
photons in the assigned wavelength range 



7(X)AX = S[l/ M (X)](AX/<$» 



(1) 



per unit area and per unit time. This approximate 
spectral distribution is expressed in the scale of 
wavelengths because it is in this scale that the 
statistical disl ribution of successive steps of degrada- 
tion is more nearly a constant. 

Equation (1) is subject to three main limitations: 
(a) It does not bold at photon energies > 500 kev. 
where the distribution of degradation steps varies 
rapidly from step to step; (b) it represents the 
spectrum for a "steady state" of degration, which 
does not hold for photons that have been scattered 
only a very few times because their spectral distri- 
bution depends on the source wavelength ("tran- 
sient" effect) (see footnote 4); (c) it does not take 
into account outright photon absorption, by photo- 
electric effect. The effect of absorption can easily 
be included in (1) as a progressive decay of intensity 
as X increases, provided the probability of absorption 
at each step is small. 

Equation (1) is contained as a special case in the 
solution of degradation problems by successive ap- 
proximations (see footnote 4). In this solution the 
degradation spectrum is resolved into steady state 
and transient components. A numerical calculation 
of the first few terms of this solution for the X-ray 
degradation would be possible, but it was found 
more convenient to proceed as in the analogous 
electron degradation problem. 7 One starts from an 
integral form of the degradation equation, namely 

£d\'I(\')Li a * 9 (\')+ f V+a ifc(X / ,X // )dX // ]=S (2) 

which represents a sort of "balance sheet" of photons 
in the entire wavelength range lrom X to X. The 
source injects S photons per unit volume and unit 

time at X . Of these, rfX / /(X / )Mabs(X / ) are lost 

through outright absorption, and the remainder by 
Compton scattering with final wavelength X">X, 

7 L. V. Spencer and U. Fano, Phys. Rev. 93, 1172 (1954). 
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where &(X',X")^X" is the probability per unit path- 
length that a photon of wavelength X' be scattered 
to a final wavelength between X" and X"+g?X". 
Note that the wavelengths are expressed in Compton 
units of h/mc=2AXlO~ 10 cm. 

Equation (2) should now be transformed so that 
it displays an approximate solution similar to (1) 
and so that the necessary corrections can be easily 
estimated and then calculated either by successive 
approximation or by numerical solution of the new 
equation. To minimize the difficulties (a) and (b) 
indicated above, we restrict our discussion to wave- 
lengths X>X +2; thus no photon from the source can 
reach X in a single Compton scattering, because 
&(X',X") = for X">X'+2. As a first step we trans- 
form the Compton scattering term in (2), by setting 
X" = X'-H and X' = X— rj, and find 

rw^x') r +2 yt(x , ,x ,/ )c/x ,/ 

Jx J\ 

Jo J v 

= f# (*dnI(k-r,)H\-r,,\-r, + t). (3) 

Jo Jo 

Because the procedure that is being developed is 
expected to apply when / and k vary slowly with 77, 
it is convenient to represent (3) in a form suited to 
expansion in powers of 77. This is done by an 
operator symbol based on the power expansion 
formula 

j{\- r,) =/(X) - „ £ /(\)+i f ^ /(X) 



d\ 



dX 2 



15S , 



+ ..=« ^j{\). (4) 



The integrand on the right of (3) can therefore be 
indicated in the form 

_ d_ 
I(\~r ] )k(\-v^-v+^=e Vd *I(\)k(\,\+$. (5) 

We consider in particular the value of (3) for 77 set 
equal to zero in the exponent of (5). This value 
is 



/(x)J q 2 ^ 6fc(x, 



X+f)=/(X)M 1 (X). 



(6) 



The function Mi(X) defined here is the first moment 
of the wavelength shift and may be expressed as 
Mscatt(£), which coincides with the denominator of (1). 
We add the right side of (6) to the expression (3) 
and subtract its left side from (3), utilize (5) and 
combine the result with the degradation equation 
(2), which becomes 

/(X)M 1 (X)+ fWj(X'Ws(V) 

Jx 

+ f # PdXr^--l)/(X)fc(X, \+Q = S. (7) 

Jo Jo 



The solution (1) is thereby obtained if the integrals 
in this equation can be disregarded by assuming 
Mabs ~ 0, and 7/d/dX ^ in the exponent. 

Consider now the absorption, represented by the 
first integral in (7), taking into account that 1M X 
should be effectively constant in the absence of 
absorption. The absorption term is then written 

I V M abs(X / )/(X / )= ( V^xo^xwcx')], (8) 

J \ J x o 

showing that IM X is effectively reduced, as X in- 
creases, by a fraction 

0(X)AX=^AX=-^W^ (9) 

MlQ<) MscattW (0 

in every small interval AX. Notice that this fraction 
has the dimensions of a pure number. This absorp- 
tion effect, taken by itself, would cause IM X to 

decay, as X increases, according to exp — (X') d\' • 

Accordingly we replace the spectral distribution 
I(\) with a new unknown w(\) ~ 1 defined by the 
expression 



/(X) = 



S - f CM**' 

p JXo 



"M 1 (X) t 



w(X). 



(10) 



When this expression for 7(X) is substituted in (7), 
the source strength S drops out as a common factor. 
In the first term we use the identity 

px' r\> 

e Jxo =1 _ d\'l3(\')e Jx ° , (11) 

JXo 

so that the second term from this expression may be 
combined with the absorption term in (7). In the 

second integral of (7) the term exp — /3(\')d\' 

may be factored out utilizing the operator formula 



"ax 



-f 0(X')dX' -f /3(X')dX' -^[^-^(x)] 

f JX -_,, JXo p LdX J. 



; (12) 



after factoring, (11) is applied again. Equation 
(7) becomes now 

rx -( /s(x")rfx" r rk rz 

w(\)+\ dk'p(\')e Jx ° < w (X')- \ di \ dri 

JXo ^ Jo Jo 



'[sh'wiAtx.x+a 



M,(X) 



w(X) 



J Jo Jo 



dr\ 






x+|) 

(X) 



w(X)=l. 



(13) 



This equation is fully equivalent to the initial 
eq (2), for X>Xo+2. Its zero-order approximation 
solution is w(X) = l, obtained by disregarding all the 
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integrals. To obtain a better approximation, the 
integrals may be evaluated by whatever procedure 
appears expedient. We proceed here by expanding 
into powers of the operator d/dX and of the rate of 
absorption (3, to first order only. The integral 
over X' yields no contribution to first order, because 
of the factor /3(X') and because the factor contained 
in the curly brackets vanishes in zero-th order. 
Where the zero order approximation { }=0 is poor, 
for X<X + 2, /3(X') is very small for most materials. 
The second integral is evaluated to first order by 
writing 



/>(-- [ *- w] -o-j><te- 



m~\ 



(14) 



and w(X)»l. Integration over £ yields the first 
order estimate of the second integral 



-|_^ () J 2 m;oo (15) 

where M 2 (X)= I d£ £ 2 &(X,X + £) indicates the second 
moment of the wavelength shift £. Thus we find 



M»=i+£\m-K>-) iu > c 



dX 2 Mx(X) 



and 



/(X)= 



s 



MtiX) 



-J: 



v,» 



,)dv [l+ 



2 M t (X) 



b 1 M 2 (X) 
dX2M,(X) 

1 M 2 (X) 

2 M^xy 



(16) 



•]• 



(17) 



The function (1/2)M 2 /Mi was determined analyti- 
cally using the Klein-Nishina formula for Compton 
scattering. A graphical differentiation was used to 
obtain (d/5X)(l/2)M 2 /Mi. These data, which are 
the same for all materials, are given in figure 1. 

The absorption exponent dX'j8(X'), where only 

photoelectric absorption is included in the tabulations 
of jft(X'), is given for a few materials in figure 2. 
The computation of jS(X') utilized total attenuation 
coefficients from G. Grodstein 8 and the total 
Compton cross section from the Klein-Nishina 
formula. The integration was carried out numeri- 
cally. Figure 3 gives values of Mi" 1 (X). The solid 
curve is a dimensionless plot of the quantity 
[(Mcompton/MThompson)(?)]~ 1 , which is independent of 
atomic number. For a particular material, the 
nomograph on the left permits the conversion to 
units of g/cm 2 . If the source yields S photons 
g -1 sec -1 , S/Mi(X) is in units cm~ 2 sec _1 (h/mc) _1 . 

With the data from figures 1 to 3 the zero-th order 
expression may be calculated, and the importance 
of the first order correction may be evaluated 

8 G. W. Grodstein, NBS Circular 583 (1957). 
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Figure 1. Graph of corrective parameters in eq 17. 
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FlGURE 2. 
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Figure 3. Spectral distribution of photon track length. 
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Figure 4. Degradation spectra for LiH and H 2 0. 
, Approximation, eq (17); O.A> exact numerical integration. 

readily. Tliis has been done for two eases, namely 
gamma rays with initial wavelength Xo=l in LiH 
and H 2 0, the source strength S being unity. The 
results are the curves in figure 4. The points in the 
figure show the results of direct numerical integration. 
The agreement is well within the expected accuracy 
of the approximation. The departures at increasing 
wavelength reflect the increasing strength of photo- 
absorption, which leads eventually to a breakdown 
of the validity of the assumption that the absorption 
changes little over a wavelength range (£)~1. 

Of general interest is the absorption "cut-off" of 
the spectra. This point, which may be characterized 

by a parameter \ z satisfying the relation I z /3(\') 
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Figure 5. Absorption "cut-off' of the spectral distribution. 

dy = l. It is the wavelength or energy at which 
the absorption reduces sharply the spectral distribu- 
tion. Figure 5 gives a curve of \ z versus the atomic 
number Z. Again, the values for jS(X') take into 
account only photoelectric absorption. 
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